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We study the effect of continuous quantum error correction in the case where each qubit in 
a codeword is subject to a general Hamiltonian interaction with an independent bath. We first 
consider the scheme in the case of a trivial single-qubit code, which provides useful insights into the 
workings of continuous error correction and the difference between Markovian and non-Markovian 
decoherence. We then study the model of a bit-flip code with each qubit coupled to an independent 
bath qubit and subject to continuous correction, and find its solution. We show that for sufficiently 
large error-correction rates, the encoded state approximately follows an evolution of the type of a 
single decohering qubit, but with an effectively decreased coupling constant. The factor by which the 
coupling constant is decreased scales quadratically with the error-correction rate. This is compared 
to the case of Markovian noise, where the decoherence rate is effectively decreased by a factor 
which scales only linearly with the rate of error correction. The quadratic enhancement depends on 
the existence of a Zeno regime in the Hamiltonian evolution which is absent in purely Markovian 
dynamics. We analyze the range of validity of this result and identify two relevant time scales. 
Finally, we extend the result to more general codes and argue that the performance of contirmous 
error correction will exhibit the same qualitative characteristics. 



I. INTRODUCTION 



Reliable information processing requires the ability to store and manipulate information with practically negligible 
loss. Information carriers, however, constantly interact with their surroundings, which poses the risk of information 
being irreversibly dissipated. This problem is of particular significance in the case of quantum information, due 
to the inherent fragility of quantum superpositions in the presence of external interactions. Such interactions can 
quickly lead to entanglement between the system of interest and its environment, effectively resulting in the loss of 
information. This process, known as decoherence, is a major obstacle in the construction of large-scale quantum 
information devices, since as quantum systems grow in size, they also become increasingly difficult to isolate from 
their environment. 

Even though decoherence may seem to be a fundamental difficulty, the development of the theory of quantum fault 
tolerance [ll, |3, S 0j [E| has shown that it is possible in principle to implement reliable quantum information processing 
with systems of any size. As long as the error rate per information unit per time step is kept below a certain threshold, 
quantum information can be processed with an arbitrarily small error. This result is based on the idea of quantum 
error correction [1, d, 0], where the quantum state of a single information unit, say a qubit, is encoded in the state 
of a larger number of qubits. The encoding is such that if a single qubit in the code undergoes an error, the original 
state can be recovered by applying an appropriate measurement on the codeword followed by a correcting operation. 
The success of this scheme depends on the assumption that individual qubits undergo independent errors with small 
probability, and thus that errors on multiple qubits have probabilities of higher order. This technique can be extended 
to multi-qubit errors by constructing more complicated codes or by concatenation Q. 
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A. Continuous quantum error correction 



In general, error probabilities increase with time. No matter how complicated a code or how many levels of 
concatenation are involved, the probability of uncorrectable errors is never truly zero, and if the system is exposed to 
noise for a sufficiently long time the weight of uncorrectable errors can accumulate. To combat this, error correction 
must be applied repeatedly and sufficiently often. If one assumes that the time for an error-correcting operation is small 
compared to other relevant time scales of the system, error-correcting operations can be considered instantaneous. 
Then the scenario of repeated error correction leads to a discrete evolution which often may be difficult to describe. 
To study the evolution of a system in the limit of frequently applied instantaneous error correction, Paz and Zurek 
proposed to describe error correction as a continuous quantum jump process @. In this model, the infinitesimal 
error-correcting transformation that the density matrix of the encoded system undergoes during a time step dt is 

p ^ {1 - Kdt)p + Kdt^{p), (1) 

where $(p) is the completely positive trace-preserving (CPTP) map describing a full error-correcting operation, and 
K is the error-correction rate. The full error-correcting operation ^(p) consists of a syndrome detection, followed (if 
necessary) by a unitary correction operation conditioned on the syndrome. 

Consider, for example, the three-qubit bit-flip code whose purpose is to protect an unknown qubit state from bit-flip 
(Pauli X) errors. The code space is spanned by |0) — |000) and |1) = and the stabilizer generators are ZZI 

and IZZ. Here by X, Y, Z and / we denote the usual Pauli operators and the identity, respectively, and a string of 
three operators represents the tensor product of operators on each of the three qubits. The standard error-correction 
procedure involves a measurement of the stabilizer generators, which projects the state onto one of the subspaces 
spanned by |000) and |100) and |011), |010) and |101), or |001) and |110); the outcome of these measurements 

is the error syndrome. Assuming that the probability for two- or three-qubit errors is negligible, then with high 
probability the result of this measurement is cither the original state with no errors, or with a single X error on the 
first, the second, or the third qubit. Depending on the outcome, one then applies an X gate to the erroneous qubit 
and transforms the state back to the original one. The CPTP map ^{p) for this code can be written explicitly as 

$(p) = (|000)(000| + |lll>(lll|)p(|000)(000|-l- 1111>(111|) 
+ (|000)(100| + |lll>(011|)p(|100)(000| + |011>(111|) 

+ (|000>(010| + |lll)(101|)p(|010>(000| + |101>(lll|) ^ ' 

-H(|000>(001| + |lll>(110|)p(|001>(000| + |110>(111|) 

The quantum-jump process ([T]) can be viewed as a smoothed version of the discrete scenario of repeated error 
correction, in which instantaneous full error-correcting operations are applied at random times with rate k. It can 
also be looked upon as arising from a continuous sequence of infinitesimal CPTP maps of the type |T]). In practice, 
such a weak map is never truly infinitesimal, but rather has the form 

p^(l-e)p + e$(p), (3) 

where e ^ 1 is a small but finite parameter, and the weak operation takes a small but nonzero time Tc. For times t 
much greater than Tc {tc <^ t), the weak error-correcting map ((3]) is well approximated by the infinitesimal form p]). 
where the rate of error correction is 

K = e/tc. (4) 

A weak map of the form ([3]) could be implemented, for example, by a weak coupling between the system and an ancilla 
via an appropriate Hamiltonian, followed by discarding the ancilla. A closely related scenario, where the ancilla is 
continuously cooled in order to reset it to its initial state, was studied in [loj . 

Another way of implementing the weak map is via weak measurements followed by weak unitaries dependent on 
the outcome. The corresponding weak measurements, however, are not weak versions of the strong measurements for 
syndrome detection; they are in a different basis [11]. They can be regarded as weak versions of a different set of 
strong measurements which, when followed by an appropriate unitary, yield the same map ^(p) on average. Thus, 
the workings of continuous error correction, when it is driven by weak measurements, does not translate directly into 
the error syndrome detection and correction of the standard paradigm. In this sense, the continuous approach can be 
regarded as a different paradigm for error correction — one based on weak measurements and weak unitary operations. 
The idea of using continuous weak measurements and unitary operations for error correction has been explored in the 
context of different heuristic schemes [l^ ITsj , some of which are based on a direct "continuization" of the syndrome 
measurements. In this paper we consider continuous error correction of the type given by Eq. P]). 
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B. Markovian decoherence 



So far, continuous quantum error correction has been studied only for Markovian error models. The Markovian 
approximation describes situations where the bath-correlation times are much shorter than any characteristic time 
scale of the system [13] • In this limit, the dynamics can be described by a semi-group master equation in the Lindblad 
form 

^ = L{p) EE p] + i^A, (2L,pLt - l]LjP - pL]L,). (5) 

3 

Here H is the system Hamiltonian and the {Lj} are suitably normalized Lindblad operators describing different error 
channels with decoherence rates \j. For example, the Liouvillian 

L{p)=Y,\,{X,pX,-p), (6) 

3 

where Xj denotes a local bit-flip operator acting on the j'-th qubit, describes independent Markovian bit-flip errors. 

For a system undergoing Markovian decoherence and error correction of the type ^ , the evolution is given by the 
equation 

g^L(p) + ^r(p), (7) 

where r{p) = — p- In 9], Paz and Zurek showed that if the set of errors {Lj} are correctable by the code, in the 
limit of infinite error-correction rate (strong error-correcting operations applied continuously often) the state of the 
system freezes and is protected from errors at all times. The effect of freezing can be understood by noticing that the 
transformation arising from decoherence during a short time step At, is 

p + L{p)At+ 0{Ae), (8) 

i.e., the weight of correctable errors emerging during this time interval is proportional to At, whereas uncorrectable 
errors (e.g. multi-qubit bit flips in the case of the thrcc-qubit bit-flip code) are of order 0{At^). Thus, if errors are 
constantly corrected, in the limit At uncorrectable errors cannot accumulate, and the evolution stops. 



C. The Zeno effect. Error correction versus error prevention 

The effect of "freezing" in continuous error correction strongly resembles the quantum Zeno effect (l6j . in which 
frequent measurements slow down the evolution of a system, freezing the state in the limit where they are applied 
continuously. The Zeno effect arises when the system and its environment are initially decoupled and they undergo 
a Hamiltonian-driven evolution, which leads to a quadratic change with time of the state during the initial moments 
[l3| (the so called Zeno regime). Let the initial state of the system plus the bath be psb(O) = |0)(0|5 (g) pb(0). For 
small times, the fidelity of the system's density matrix with the initial state a{t) = Tr {(|0) {0\s (8) Ib) Pssit)} can be 
approximated as 

a{t) ^1-Ct^ + 0{t^). (9) 

In terms of the Hamiltonian Hs b acting on the entire system, the coefficient C is 

C = Tr {i/ls (|0)(0|5 ® Pb(0))} - Tv{HsB (|0)(0|s ® Ib) Hsb (|0)(0|s (g> pb(0))} . (10) 

According to Eq. if after a short time step At the system is measured in an orthogonal basis which includes the 
initial state |0), the probability to find the system in a state other than the initial state is of order 0{At^). Thus if 
the state is continuously measured {At 0), this prevents the system from evolving. 

It has been proposed to utilize the quantum Zeno effect in schemes for error prevention [isl . [l^ [20l | , in which an 
unknown encoded state is prevented from errors simply by frequent measurements which keep it inside the code space. 
The approach is similar to error correction in that the errors for which the code is designed send a codeword to a space 
orthogonal to the code space. The difference is that different errors need not be distinguishable, since the procedure 
does not involve correction of errors, but their prevention. In [20| it was shown that with this approach it is possible 
to use codes of smaller redundancy than those needed for error correction and a four-qubit encoding of a qubit was 
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proposed, which is capable of preventing arbitrary independent errors arising from Hamiltonian interactions. The 
possibility of this approach implicitly assumes the existence of a Zeno regime, and fails if we assume Markovian 
decoherence for all times. This is because the probability of errors emerging during a time step dt in a Markovian 
model is proportional to dt (rather than dt^), and hence errors will accumulate with time if not corrected. 

From the above observations we see that error correction is capable of achieving results in noise regimes where error 
prevention fails. Of course, this advantage is at the expense of a more complicated procedure — in addition to the 
measurements used in error prevention, error correction involves unitary correction operations, and in general requires 
codes with higher redundancy. At the same time, we see that in the Zeno regime it is possible to reduce decoherence 
using weaker resources than those needed in the case of Markovian noise. This suggests that in this regime error 
correction may exhibit higher performance than it does for Markovian decoherence. 

D. Non-Mar kovian decoherence 

Markovian decoherence is an approximation valid for times much larger than the memory of the environment. 
In many situations of practical significance, however, the memory of the environment cannot be neglected and the 
evolution is highly non-Markovian [3, HH [11, . Furthermore, no evolution is strictly Markovian, and for a system 
initially decoupled from its environment a Zeno regime is always present, short though it may be [l7| . If the time 
resolution of error-correcting operations is high enough so that they "see" the Zeno regime, this could give rise to 
different behavior. 

The existence of a Zeno regime is not the only interesting feature of non-Markovian decoherence. The mechanism 
by which errors accumulate in a general Hamiltonian interaction with the environment may differ significantly from 
the Markovian case, since the system may develop nontrivial correlations with the environment. For example, imagine 
that some time after the initial encoding of a system, a strong error-correcting operation is applied. This brings the 
state inside the code space, but the state contains a nonzero portion of errors non-distinguishable by the code. Thus 
the new state is mixed and is generally correlated with the environment. A subsequent error-correcting operation 
can only aim at correcting errors arising after this point, since the errors already present inside the code space are 
in principle uncorrectable. Subsequent errors on the density matrix, however, may not be completely positive due to 
the correlations with the environment. 

Nevertheless, it follows from a result in (2^ that an error-correction procedure which is capable of correcting a 
certain class of completely positive (CP) maps, can also correct any linear noise map whose operator elements can 
be expressed as linear combinations of the operator elements in a correctable CP map. This implies, in particular, 
that an error-correction procedure that can correct arbitrary single-qubit CP maps can correct arbitrary single-qubit 
linear maps. The effects of system-environment correlations in non-Markovian error models have also been studied 
from the perspective of fault tolerance, and it has been shown that the threshold theorem can be extended to various 
types of non-Markovian noise [1^, [2^ [23| . 

Another important difference from the Markovian case is that error correction and the effective noise on the 
reduced density matrix of the system cannot be treated as independent processes. One could derive an equation for 
the effective evolution of the system alone subject to interaction with the environment, like the Nakajima-Zwanzig 
[28l [IPl or the time-convolutionless (TCL) [sol, 'ST'l master equations, but the generator of transformations at a given 
moment in general will depend (implicitly or explicitly) on the entire history up to this moment. Therefore, adding 
error correction can nontrivially affect the effective error model. This means that in studying the performance of 
continuous error correction one either has to derive an equation for the effective evolution of the encoded system, 
taking into account error correction from the very beginning, or one has to look at the evolution of the entire system — 
including the bath — where the error generator and the generator of error correction can be considered independent. 
In the latter case, for sufficiently small Tc, the evolution of the entire system including the bath can be described by 

^ = -^[H,p]+nT{p), (11) 

where p is the density matrix of the system plus bath, H is the total Hamiltonian, and the error-correction generator 
r acts locally on the encoded system. In this paper, we take this approach for a sufficiently simple bath model which 
allows us to find a solution for the evolution of the entire system. 

E. Plan of this paper 

The rest of the paper is organized as follows. To develop understanding of the workings of continuous error 
correction, in Sec. II we look at a simple example: an error-correction code consisting of only one qubit which aims at 
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protecting a known state. We discuss the difference in performance for Markovian and non-Markovian decoherence, 
and argue the impUcations it has for the case of multi-qubit codes. In Sec. Ill, we study the three-qubit bit-flip code. 
We first review the performance of continuous error correction in the case of Markovian bit-flip decoherence, which 
was first studied in We then consider a non-Markovian model, where each qubit in the code is coupled to an 
independent bath qubit. This model is sufficiently simple so that we can solve for its evolution analytically. In the 
limit of large error-correction rates, the effective evolution approaches the evolution of a single qubit without error 
correction, but the coupling strength is now decreased by a factor which scales quadratically with the error-correction 
rate. This is opposed to the case of Markovian decoherence, where the same factor scales linearly with the rate of 
error-correction. In Sec. IV, we show that the quadratic enhancement in the performance over the case of Markovian 
noise can be attributed to the presence of a Zeno regime and argue that for general stabilizer codes and independent 
errors, the performance of continuous error correction would exhibit the same qualitative characteristics. In Sec. V, 
we conclude. 



II. THE SINGLE-QUBIT CODE 



Consider the problem of protecting a qubit in state |0) from bit- flip errors. This problem can be regarded as a 
trivial example of a stabilizer code, where the code space is spanned by |0) and its stabilizer is Z. Let us consider the 
Markovian bit-fiip model first. The evolution of the state subject to bit-flip errors and error correction is described 
by Eq. ^ with 

L{p) - X{XpX - p), (12) 

and 

r(p) = |o)(o|p|o)(o| + |o)(i|p|i)(o|-p. (13) 

If the state lies on the z-axis of the Bloch sphere, it will never leave it, since both the noise generator and the 
error-correction generator (fT5|) keep it on the axis. We will take the qubit to be initially in the desired state |0), and 
therefore at any later moment it will have the form p{t) = a{t)\0){0\ -I- (1 — a(t))|l)(l|, a{t) G [0, 1]. The coefficient 
a(t) has the interpretation of a fidelity with the trivial code space spanned by |0). For an infinitesimal time step dt, 
the effect of the noise is to decrease a{t) by the amount X{2a{t) — l)dt and that of the correcting operation is to 
increase it by k(1 — a(t))dt. The net evolution is then described by 

^^^~{K + 2X)a{t) + {K + X). (14) 
dt 

The solution is 

a(<) = (l-af)e-(''+2^)* + «f, (15) 

where 

af = l--^, (16) 
2 + r 

and r = k/A is the ratio between the rate of error correction and the rate of decoherence. We see that the fidelity 
decays, but it is confined above its asymptotic value a^, which can be made arbitrarily close to 1 for a sufficiently 
large r. 

Now let us consider a non-Markovian error model. We choose the simple scenario where the system is coupled to 
a single bath qubit via the Hamiltonian 

H = -fX(g>X, (17) 

where 7 is the coupling strength. This can be a good approximation for situations in which the coupling to a single 
spin from the bath dominates over other interactions [23| . 

We will assume that the bath qubit is initially in the maximally mixed state, which can be thought of as an 
equilibrium state at high temperature. From Eq. (jlip one can verify that if the system is initially in the state |0), 
the state of the system plus the bath at any moment will have the form 

Pit) = (a(i)|0)(0| + (1 - a(i))|l)(l|) ® ^ - P{t)Y ® ^. (18) 
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In the tensor product, the first operator belongs to the Hilbert space of the system and the second to the Hilbert 
space of the bath. We have a{t) £ [0, 1], and \f3{t)\ < y^a{t){l — a{t)),(3{t) € R. The reduced density matrix of the 
system has the same form as the one for the Markovian case. The traceless term proportional to f3{t) can be thought 
of as a "hidden" part, which nevertheless plays an important role in the error-creation process, since errors can be 
thought of as being transferred to the "visible" part from the "hidden" part (and vice versa). This can be seen from 
the fact that during an infinitesimal time step dt, the Hamiltonian changes the parameters a and (3 as follows: 



(3-^ 13+ {2a- l)jdt. 

The effect of an infinitesimal error-correcting operation is 

a ^ a -|- (1 — a)Kdt, 
(3^ (3- (3Kdt. 

Note that the hidden part is also being acted upon. Putting it all together, we get the system of equations 

da{t) 



(19) 



(20) 



dt 

d(3{t) 
dt 



nil 



a{t))-2^m, 
7(2a - 1) - K/3(i). 



The solution for the fidelity a{t) is 



a{t) = 



27^ 



47^ 



K7 



47^ 



■ sin 27^ - 



27^ 



4^2 



■ cos 2^t 



We see that as time increases, the fidelity stabilizes at the value 

2 + i?2 



„NM 



i?2 



= 1 - 



(21) 



(22) 



(23) 



where R = k/j is the ratio between the error-correction rate and the coupling strength. In Fig. 1 we have plotted 
the fidelity as a function of the dimensionless parameter jt for three different values of R. For error-correction rates 
comparable to the coupling strength {R — 1), the fidelity undergoes a few partial recurrences before it stabilizes close 
to a^^. For larger i? = 2, however, the oscillations are already heavily damped and for i? = 5 the fidelity seems 
confined above a^^. As R increases, the evolution becomes closer to a decay like the one in the Markovian case. 



a 



1 

. 9 



. 7 
. 6 

. 4 





\\ 








^ \ 








\ . 

\ ^. 

\ ^ _ . 








\ 

\ 

\ 

\ 








\ 














— - R=l 
-- R=2 
— R=5 


1 


2 


3 



yt 



Figure 1: (Color online) Fidelity of the single-qubit code with continuous bit-flip errors and correction, as a function of 
dimensionless time ^t, for three different values of the ratio R = n/^. 



A remarkable difference, however, is that the asymptotic weight outside the code space (1 — a^^) decreases with k 
as 1/k^, whereas in the Markovian case the same quantity decreases as 1/k. The asymptotic value can be obtained 
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as an equilibrium point at which the infinitesimal weight fiowing out of the code space during a time step dt is equal 
to the weight flowing into it. The latter corresponds to vanishing right-hand sides in Eqs. (|14p and (HI]). In Sec. IV, 
we will show that the difference in the equilibrium code-space fidelity for the two different types of decoherence arises 
from the difference in the corresponding evolutions during initial times. 

For multi-qubit codes, error correction cannot preserve a high fidelity with the initial codeword for all times, because 
there will be multi-qubit errors that can lead to errors within the code space itself. But it is natural to expect that the 
code-space fidelity can be kept above a certain value, since the effect of the error-correcting map ([1]) is to oppose its 
decrease. If similarly to the single-qubit code there is a quadratic difference in the code-space fidelity for the cases of 
Markovian and non-Markovian decoherence, this could lead to a different performance of the error-correction scheme 
with respect to the rate of accumulation of uncorrectable errors inside the code space. This is because multi-qubit 
errors that can lead to transformations entirely within the code space during a time step dt are of order 0{dt^). This 
means that if the state is kept constantly inside the code space (as in the limit of an infinite error-correction rate), 
uncorrectable errors will never develop. But if there is a finite nonzero portion of correctable errors, by the error 
mechanism it will give rise to errors not distinguishable or misinterpreted by the code. Therefore, the weight outside 
the code space can be thought of as responsible for the accumulation of uncorrectable errors, and consequently a 
difference in its magnitude may lead to a difference in the overall performance. In the following sections we will see 
that this is indeed the case. 



III. THE THREE-QUBIT BIT-FLIP CODE 
A. A Markovian error model 



Even though the three-qubit bit-fiip code can correct only bit-fiip errors, it captures most of the important charac- 
teristics of nontrivial stabilizer codes. Before we look at a non-Markovian model, we will review the Markovian case 
which was studied in [§]. Let the system decohere through identical independent bit-fiip channels, i.e., L{p) is of the 
form ([6]) with Ai = A2 = A3 = A. Then one can verify that the density matrix at any moment can be written as 



where 



p{t) = a(<)p(0) + b{t)pi + c{t)p2 + d{t)p3, (24) 
pi = ^{X,p{0)X, + X2p{0)X2 + X3piO)X3), 

P2 - ^{X,X2P{0)X,X2 + X2X3piO)X2X3 + X,X3piO)X,X3), (25) 
P3 = XiX2X3p{0)XiX2X3, 

are equally-weighted mixtures of single-qubit, two-qubit and three-qubit errors on the original state. 

The effect of decoherence for a single time step dt is equivalent to the following transformation of the coefficients 
in Eq. ((M)) : 

a — > a — 3aXdt + bXdt, 
b^b + iaXdt - 3bXdt + 2cXdt, 

c^c + 2bXdt - 3cXdt + MXdt, ^ ^ 

d ~f d + cXdt — SdXdt. 



If the system is initially inside the code space, combining Eq. ((26)) with the effect of the weak error-correcting map 
p ^ (1 — Kdt)p + Kdt^{p), where $(p) is given in Eq. yields the following system of first-order linear differential 
equations for the evolution of the system subject to decoherence plus error correction: 

^ = -3Aa(t) + (A + K)&(t), 
^ = 3Xa{t) - (3A + K)bit) + 2Xc{t), 

dm ^^^^ 
^ = 2Xb{t) - (3A + K)c(i) + 3Xd{t), 

dd{t) 



dt 



(A -I- K)c{t) - iXd{t). 
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The exact solution has been found in Here we just note that for the initial conditions a(0) = 1,6(0) = c(0) = 
d{0) — 0, the exact solution for the weight outside the code space is 

6(t)+c(t) = -^(l-e-(4+'-)^*), (28) 
4 + 7- 

where r — k/X. We see that similarly to what we obtained for the trivial code in the previous section, the weight 
outside the code space quickly decays to its asymptotic value which scales as l/r. But note that here the 
asymptotic value is roughly three times greater than that for the single-qubit model. This corresponds to the fact 
that there are three single-qubit channels. More precisely, it can be verified that if for a given k the uncorrected 
weight by the single-qubit scheme is small, then the uncorrected weight by a multi-qubit code using the same k and 
the same kind of decoherence for each qubit scales approximately linearly with the number of qubits [ll| . Similarly, 
the ratio r required to preserve a given overlap with the code space scales linearly with the number of qubits in the 
code. 

The most important difference from the single-qubit model is that in this model there are uncorrectable errors that 
cause a decay of the state's fidelity inside the code space. Due to the finiteness of the resources employed by our 
scheme, there always remains a nonzero portion of the state outside the code space, which gives rise to uncorrectable 
three-qubit errors. To understand how the state decays inside the code space, we ignore the terms of the order of the 
weight outside the code space in the exact solution. We obtain: 

«(0«^^ «l-d(t), (29) 



b{t) w c{t) w 0. (30) 

Comparing this solution to the expression for the fidelity of a single decaying qubit without error correction — which 
can be seen from Eq. p3|) for k — — we see that the encoded qubit decays roughly as if subject to bit-flip decoherence 
with rate 6A/r. Therefore, for large r this error-correction scheme can reduce the rate of decoherence approximately 
r/6 times. In the limit r —^ oo, it leads to perfect protection of the state for all times. 



B. A non-Markovian error model 



We consider a model where each qubit independently undergoes the same kind of non-Markovian decoherence as 
the one we studied for the single-qubit code. Here the system we look at consists of six qubits - three for the codeword 
and three for the environment. We assume that all system qubits are coupled to their corresponding environment 
qubits with the same coupling strength, i.e., the Hamiltonian is 

3 

H^^J2^f<E>Xf, (31) 

i=l 

where the operators act on the system qubits and act on the corresponding bath qubits. The subscripts label 
the particular qubit on which they act. Obviously, the types of effective single-qubit errors on the density matrix 
of the system that can result from this Hamiltonian at any time, whether they are CP or not, will have operator 
elements which are linear combinations of / and X^ , i.e., they are correctable by the procedure according to [24j |. 
Considering the forms of the Hamiltonian (PTjl and the error-correcting map , one can see that the density matrix 
of the entire system at any moment is a linear combination of terms of the following type: 

vl+P -y-m+q -yn+r 

Qi^r.,pqr = X[X^XSp{Q)XfXlX-, (g> ® ® (32) 

Here the first term in the tensor product refers to the Hilbert space of the system, and the following three refer to 
the Hilbert spaces of the bath qubits that couple to the first, second and third qubits from the code, respectively. 
The powers Z,m,n,p, g, r take values and 1 in all possible combinations, and X^ = X, X^ = X^ ~ I. Note that 
Qimn,pqr should not be mistaken for the components of the density matrix in the computational basis. Collecting these 
together, we can write the density matrix in the form 

pit) — ^ ^ ( *) ^ (i)^''''''' Clmn,pqr{t) ^ Qlmn.pqn (33) 
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where the coefBcients Cimn,pqr{t) are real. The coefficient Cooo,ooo is less than or equal to the codeword fidelity (with 
equality when p{0) — |0)(0| or p{0) = Since the scheme is intended to protect an unknown codeword, we are 

interested in its worst-case performance; we will therefore use Cooo,ooo a-s a lower bound on the codeword fidelity. 

Using the symmetry with respect to permutations of the different system-bath pairs of qubits and the Hermiticity 
of the density matrix, we can reduce the description of the evolution to a system of equations for only 13 of the 64 
coefficients. (In fact, 12 coefficients are sufficient if we invoke the normalization condition Trp = 1, but we have found 
it more convenient to work with 13.) The equations are linear, and we write them as a single 13-dimensional vector 
equation: 



dt 



Cooo.ooo 
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Cm, 111 



(34) 



where R — k/^. Each nonzero component in this matrix represents an allowed transition process for the quantum 
states; these transitions can be driven either by the decoherence process or the continuous error-correction process. 
We plot these allowed transitions in Fig. 2. 

We can use the symmetries of the process to recover the 64 coefficients of the full state. Each of the 13 coefficients 
represents a set of coefficients having the same number of Is on the left and the same number of Is on the right, 
as well as the same number of places which have 1 on both sides. All such coefficients are equal at all times. For 
example, the coefficient Ciio,oii is equal to all coefficients with two Is on the left, two Is on the right and exactly one 
place with 1 on both sides; there are exactly six such coefficients: 



Ci 



10,011 



= Ci 



10,101 



= C 



101,011 



= C 



101,110 



-c 



011,110 



= c 



011,101- 



In determining the transfer rate from one coefficient to another in Fig. 2, one has to take into account the number 
of different coefficients of the first type which can make a transition to a coefficient of the second type of order dt 
according to Eq. (fTTj) . The sign of the flow is determined from the phases in front of the coefficients in Eq. p3|) . 
The eigenvalues of the matrix in Eq. ([M]) up to the first two lowest orders in 1/k are presented in Table I. 



Table I: Eigenvalues of the matrix 



Eigenvalues 
Ao = 

Al,2 = —K. 

A3, 4 = — K ± i2j 

A5,6 = — K ± i47 

A7,8 = -K± ijVTS + 3)7 + 0(1/ k) 
A9,io = -K ± ijVTs ~ 3)7 + 0{1/k) 
All, 12 = ±i(24/fl^)7 ^ (144/fl^)7 + 0{1/k*) 



Obviously all eigenvalues except the first one and the last two describe fast decays with rates ^ k. They correspond 
to terms in the solution which will vanish quickly after the beginning of the evolution. The eigenvalue corresponds 
to the asymptotic {t 00) solution, since all other terms will eventually decay. The last two eigenvalues are those 
that play the main role in the evolution on a time scale t ^ ^. We see that on such a time scale, the solution will 
contain an oscillation with an angular frequency approximately equal to (24/i?^)7 which is damped by a decay factor 
with a rate of approximately (144/i?'^)7. In Fig. 3 we have plotted the codeword fidehty Cooo,ooo(i) as a function of 
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Figure 2: (Color online) These are the allowed transitions between the different components of the system (I34p and their rates, 
arising from both the decoherence (bit-flip) process (with rate 7 and the continuous error-correction process (with rate k). 
Online, the transitions due to decoherence are black, and the transitions due to error correction are red. 



the dimcnsionless parameter jt for R = 100. The graph indeed represents this type of behavior, except for very short 
times after the beginning (-yt ^ 0.1), where one can see a fast but small in magnitude decay (Fig. 4). The maximum 
magnitude of this quickly decaying term obviously decreases with R, since in the limit of R ^ 00 the fidelity should 
remain constantly equal to 1. 

From the form of the eigenvalues one can see that as R increases, the frequency of the main oscillation decreases 
as 1/R^ while the rate of decay decreases faster, as 1/R^. Thus in the limit R ^ 00, the evolution approaches 
an oscillation with an angular frequency (24/i?^)7. (We formulate this statement more rigorously below.) This is 
the same type of evolution as that of a single qubit interacting with its environment, but the coupling constant is 
effectively reduced by a factor of i?^/12. 

While the coupling constant serves to characterize the decoherence process in this particular case, this is not valid 
in general. To handle the more general situation, we propose to use the instantaneous rate of decrease of the codeword 
fidelity Few as a measure of the effect of decoherence: 



dFcJ^ 
dt 



(35) 



(In the present case, F^w — Cooo.ooo-) This quantity does not coincide with the decoherence rate in the Markovian case 
(which can be defined naturally from the Lindblad equation) , but it is a good estimate of the rate of loss of fidelity 
and can be used for any decoherence model. From now on we will refer to it simply as an error rate, but we note 
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Figure 3: (Color online) Long-time behavior of three-qubit system with bit-flip noise and continuous error correction. The ratio 
of correction rate to decoherence rate is R = K/7 — 100. 
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Figure 4: (Color online) Short-time behavior of three-qubit system with bit-flip noise and continuous error correction. The 
ratio of correction rate to decoherence rate is R = K/7 — 100. 



that there are other possible definitions of instantaneous error rate suitable for non-Markovian decoherence, which in 
general may depend on the kind of errors they describe. Since the goal of error correction is to preserve the codeword 
fidelity, the quantity ([55]) is a useful indicator for the performance of a given scheme. Note that A(_F'cw) is a function of 
the codeword fidchty and therefore it makes sense to use it for a comparison between different cases only for identical 
values of Few For our example, the fact that the coupling constant is effectively reduced approximately E? /\2 times 
implies that the error rate for a given value of F^w is also reduced B? /\2 times. Similarly, the reduction of A by the 
factor r/6 in the Markovian case implies a reduction of A by the same factor. We see that the effective reduction of 
the error rate increases quadratically with in the non-Markovian case, whereas it increases only linearly with k in 
the Markovian case. 

Now let us rigorously derive the approximate solution to this model of non-Markovian decoherence with continuous 
error correction. Assuming that 7 ^ k (or cquivalcntly, R ^ 1), the superoperator driving the evolution of the system 
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during a time step St can be written as 

5t St St 



^cst ^ ^c^st^ ldt'e^^'^''-''^C^e^^'' + j dt' j dt" e''-^''-'"^ C^e^-^'" -''^ C^e"-'' + 



Of 
St St St 

t"-t')^ 



+ j dt'j dt"J dt"'e^«(^*-*"')/:^e^«(*"'-*")/:^e^'=(*"-*')/:^e^"*' + ... (36) 

t' t" 



We have denoted the Liouvilhan by £ = £-y + where Ck.P — K,T{p), and C-yp — —i[H, p]. 

Let 7<5< ^ 1 ^ K6t. We will derive an approximate differential equation for the evolution of p{t) by looking at 
the terms of order 6t in the change of p according to Eq. ([5^ . When k = 0, we have dp/dt — C^p, so the effect 
of on the state of the system can be seen from Eq. with k taken equal to 0. By the action of exp(£„t), 

the different terms of the density matrix transform as follows: f?ooo.ooOj ^?lll.ooo^ f?iii.iii remain unchanged, £iioo,ioo ~^ 
e Pioo,ioo + (l^e )f?ooo,oooi f?iio,iio ^ e f?iio,iio + (l~e )f?iii,iiii f?iio,ooi ~* e Piio,ooi~ (l^e )£'iii,oooi 
and all other terms are changed as ^ e^'^*g. Since nSt 1, we will ignore terms of order e"'"'*. But from Eq. ([55]) 
it can be seen that all terms except pooo.ooo, f?iii,ooo, f?ooo.iii, ^?lll.lll will get multiplied by the factor e^""^* by the 
action of exp{£f^St) in Eq. p6|) . The integrals in Eq. ((36|) also yield negligible factors, since every integral either gives 
rise to a factor of order St when the integration variable is trivially integrated, or a factor of 1/k when the variable 
participates nontrivially in the exponent. Therefore, in the above approximation these terms of the density matrix 
can be neglected, which amounts to an effective evolution entirely within the code space. According to Eq. ([M)) . the 
terms g'ooo.oooj f?iii,000j fi*!!!,!!! can couple to each other only by a triple or higher application of C-y. This means that 
if we consider the expansion up to the lowest nontrivial order in 7, we only need to look at the triple integral in Eq. 

(EH). 

Let us consider the effect of exp(CSt) on Cooo,ooo- Any change can come directly only from f?iii.ooo ^nd gooo.iii- 
The first exponent e"''"* acts on these terms as the identity. Under the action of the first operator each of these 
two terms can transform to six terms that can eventually be transformed to £iooo,ooo- They are gnofloo, ^?lOl.ooo, 
f?oii,ooo, Piii,ioOj ^?lll,olOJ f?iii,ooi: and pooo,iiOj ^?ooo,lOlI f?ooo,oii, f?ioo,iiii f?oio,iii, f?ooi,iii, with appropriate factors. 
The action of the second exponent is to multiply each of these new terms by e""'* ~* After the action of the 
second , the action of the third exponent on the relevant resultant terms will be again to multiply them by a factor 
g-K(t -t ) xhus the second and the third exponents yield a net factor of e""^* ~* \ After the second and the 
third the relevant terms that we get are gooo.ooo and f?ioo.ioo, f?oio.oiOi f?ooi.ooii each with a corresponding factor. 
Finally, the last exponent acts as the identity on £iooo,ooo and transforms each of the terms ^?loo.loo: f?oio,oio, f?ooi,ooi 
into (l-e-''(^*-*"')) 

£•000,000- Counting the number of different terms that arise at each step, and taking into account 
the factors that accompany them, we obtain: 

St St St 

Cooo,ooo ^ Cooo,ooo+ /dt'/di"/dt"'(24e-«(*"''*')-36e-'^(^*-*'))Ciii,ooo + --- 



t' t" 
24 

Cooo,ooo + Cinmo^lSt + 0{St^). (37) 



Using that Cooo,ooo + (^111,111 ~ 1, in a similar way one obtains 
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Cm, 000 — > Cm, 000 — (2Cooo,ooo — + ^i^^ )• (3^) 



For times much larger than we can write the approximate differential equations 

(iCooo.ooo 24 

t^Ciii,ooo 12 

= --^7(20000,000 - 1)- (39) 

Comparing with Eq. (jl9p . we see that the encoded qubit undergoes approximately the same type of evolution as 
that of a single qubit without error correction, but the coupling constant is effectively decreased B? /12 times. The 
solution of Eq. ((39|l yields for the codeword fidelity 

l + cos(f^7i) , , 

Cooo,ooo(i) = T^^. (40) 
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This solution is valid only with precision 0(1/ R) for times jt <C R^. This is because we ignored terms whose 
magnitudes are always of order 0{1/R) and ignored changes of order 0{'-f6t/R^) per time step 6t in the other terms. 
The latter changes could accumulate with time and become of the order of unity for times jt ^ R^, which is why the 
approximate solution is invalid for such times. In fact, if one carries out the expansion p6p to fourth order in 7, one 
obtains the approximate equations 

c^C'ooo.ooo 24 72 

— -^7'-'iii,ooo — -^71^^000,000 — tj, 

dCnifloo 12 144 

^ = --^7(,20ooo,ooo - 1) - -^7L'iii,ooo, (,41j 



which yield for the fidelity 



L^000,000lrj — 7^ • l^^j 



We see that in addition to the effective error process which is of the same type as that of a single qubit, there is an 
extra Markovian bit-flip process with rate 72^ /R^. This Markovian behavior is due to the Markovian character of 
our error-correcting procedure which, at this level of approximation, is responsible for the direct transfer of weight 
between f?ooo,ooo and giii,iii, and between pm^ooo and £iooo,iii- The exponential factor explicitly reveals the range 
of applicability of solution (HO]): with precision 0{1/R), it is valid only for times jt of up to order R^. For times 
of the order of i?^, the decay becomes significant and cannot be neglected. The exponential factor may also play 
an important role for short times of up to order R, where its contribution is bigger than that of the cosine. But in 
the latter regime the difference between the cosine and the exponent is of order 0{1/R'^), which is negligible for the 
precision that we consider. 

In general, the effective evolution that one obtains in the limit of high error-correction rate does not have to 
approach a form identical to that of a single decohering qubit. The reason we obtain such behavior here is that for 
this particular model the lowest order of uncorrectable errors that transform the state within the code space is 3, and 
thrcc-qubit errors have the form of an encoded X operation. Furthermore, the symmetry of the problem ensured an 
identical evolution of the three qubits in the code. For general stabilizer codes, the errors that a single qubit can 
undergo are not limited to bit flips only. Therefore, different combinations of single-qubit errors may lead to different 
types of lowest-order uncorrectable errors inside the code space, none of which in principle has to represent an encoded 
version of the single-qubit operations that compose it. In addition, if the noise is different for the different qubits, 
there is no unique single-qubit error model to compare to. Nevertheless, we will show that with regard to the effective 
decrease in the error-correction rate, general stabilizer codes will exhibit the same qualitative performance. 



IV. RELATION TO THE ZENO REGIME 



The effective continuous evolution pop was derived under the assumption that jSt <C 1 ^ uSt. The first inequality 
implies that St can be considered within the Zeno time scale of the system's evolution without error correction. On 
the other hand, from the relation between k and Tc in (|4]) we see that Tc <C St. Therefore, the time for implementing 
a weak error-correcting operation has to be sufficiently small so that on the Zeno time scale the error-correction 
procedure can be described approximately as a continuous Markovian process. This suggests a way of understanding 
the quadratic enhancement in the non-Markovian case based on the properties of the Zeno regime. 

Let us consider again the single-qubit code from Sec. II, but this time let the error model be any Hamiltonian-drivcn 
process. We assume that the qubit is initially in the state |0), i.e., the state of the system including the bath has 
the form p(0) = |0)(0| (8)pb(0). For times smaller than the Zeno time Stz, the evolution of the fidelity without error 
correction can be described by Eq. Q . Equation ^ naturally defines the Zeno regime in terms of a itself: 

a>az = l-C6tl. (43) 

For a single time step At <C Stz, the change in the fidelity is 

a-^ a~2VcVT^At+ 0{At^). (44) 

On the other hand, the effect of error correction during a time step At is 



a a + k(1 - a)At + 0{At^), 



(45) 
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i.e., it tends to oppose the effect of decoherence. If both processes happen simultaneously, the effect of decoherence 
will still be of the form l|44p . but the coefficient C may vary with time. This is because the presence of error-correction 
opposes the decrease of the fidelity and consequently can lead to an increase in the time for which the fidelity remains 
within the Zeno range. If this time is sufficiently long, the state of the environment could change significantly under 
the action of the Hamiltonian, thus giving rise to a different value for C in Eq. (|44l) according to Eq. (fTO|) . 

Note that the strength of the Hamiltonian puts a limit on C, and therefore this constant can vary only within a 
certain range. The equilibrium fidelity a^'^ that we obtained for the error model in Sec. II, can be thought of as the 
point at which the effects of error and error correction cancel out. For a general model, where the coefficient C may 
vary with time, this leads to a quasi-stationary equilibrium. From Eqs. (|44p and (|45p . one obtains the equilibrium 
fidelity 

a^-l-^. (46) 

In agreement with what we obtained in Sec. II, the equilibrium fidelity differs from 1 by a quantity proportional to 
1/k^. This quantity is generally quasi-stationary and can vary within a limited range. If one assumes a Markovian 
error model, for short times the fidelity changes linearly with time which leads to 1 — cx 1/k. Thus the difference 
can be attributed to the existence of a Zeno regime in the non-Markovian case. 

But what happens in the case of non-trivial codes? As we saw, there the state decays inside the code space and 
therefore can be highly correlated with the environment. Can we talk about a Zeno regime then? It turns out that 
the answer is positive. Assuming that each qubit undergoes an independent error process, then up to first order in Ai 
the Hamiltonian cannot map terms in the code space to other terms without detectable errors. (This includes both 
terms in the code space and terms from the hidden part, like £iiii,ooo in the example of the bit-flip code.) It can only 
transform terms from the code space into traceless terms from the hidden part which correspond to single-qubit errors 
(like giioo.ooo in the same example). Let |0), |1) be the two logical codewords and be an orthonormal basis that 
spans the space of all single-qubit errors. Then in the basis |0), |1), j-fAi), all the terms that can be coupled directly 
to terms inside the code space are |O)('0i|, |^/'i)(0|, |l)('0i|, From the condition of positivity of the density 

matrix, one can show that the coefficients in front of these terms are at most a{l — a) in magnitude, where a is the 
code-space fidelity. This implies that for small enough 1 — a, the change in the code-space fidelity is of the type (|44|). 
which is Zeno-like behavior. Then using only the properties of the Zeno behavior as we did above, we can conclude 
that the weight outside the code space will be kept at a quasi-stationary value of order 1/k^. Since uncorrectable 
errors enter the code space through the action of the error-correction procedure, which misinterprets some multi-qubit 
errors in the error space, the effective error rate will be limited by a factor proportional to the weight in the error 
space. That is, this will lead to an effective decrease of the error rate at least by a factor proportional to 1/k^. 

The accumulation of uncorrectable errors in the Markovian case is similar, except that in this case there is a direct 
transfer of errors between the code space and the visible part of the error space. In both cases, the error rate is 
effectively reduced by a factor which is roughly proportional to the inverse of the weight in the error space, and 
therefore the difference in the performance comes from the difference in this weight. The quasi-stationary equilibrium 
value of the code-space fidelity establishes a quasi-stationary fiow between the code space and the error space. One 
can think that this ffow effectively takes non-erroneous weight from the code space, transports it through the error 
space where it accumulates uncorrectable errors, and brings it back into the code space. Thus by minimizing the 
weight outside the code space, error correction creates a "bottleneck" which reduces the rate at which uncorrectable 
errors accumulate. 

Finally, a brief remark about the resources needed for quadratic reduction of the error rate. As pointed out above, 
two conditions are involved: one concerns the rate of error correction; the other concerns the time resolution of the 
weak error-correcting operations. Both of these quantities must be sufficiently large. There is, however, an interplay 
between the two, which involves the strength of the interaction required to implement the weak error-correcting map 
Let us imagine that the weak map is implemented by making the system interact weakly with an ancilla in 
a given state, after which the ancilla is discarded. The error-correction procedure consists of a sequence of such 
interactions, and can be thought of as a cooling process which takes away the entropy accumulated in the system as 
a result of correctable errors. If the time for which a single ancilla interacts with the system is Tc, one can verify that 
the parameter e in Eq. ([3]) would be proportional to g^r^, where g is the coupling strength between the system and 
the ancilla. From Eq. Q we then obtain that 

K (X g^Tc- (47) 

The two parameters that can be controlled are the interaction time and the interaction strength, and they determine 
the error-correction rate. Thus if g is kept constant, a decrease in the interaction time leads to a proportional 
decrease in k, which may be undesirable. In order to achieve a good working regime, one may need to adjust both 
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Tc and g. But it has to be pointed out that in some situations decreasing Tc alone can prove advantageous, if it 
leads to a time resolution revealing the non-Markovian character of an error model which was previously described 
as Markovian. The quadratic enhancement of the performance as a function of n may compensate the decrease in /t, 
thus leading to a seemingly paradoxical result: better performance with a lower error-correction rate. 

V. CONCLUSION 

In this paper we studied the performance of a particular continuous quantum error-correction scheme for non- 
Markovian errors. We analyzed the evolution of the single-qubit code and the three-qubit bit-flip code in the presence 
of continuous error correction for a simple non-Markovian bit-flip error model. This enabled us to understand the 
workings of the error-correction scheme, and the mechanism whereby uncorrectable errors accumulate. The fidelity of 
the state with the code space in both examples quickly reaches an equilibrium value, which can be made arbitrarily 
close to 1 by a sufficiently high rate of error correction. The weight of the density matrix outside the code space scales 
as 1/ K in the Markovian case, while it scales as 1/k^ in the non-Markovian case. Correspondingly, the rate at which 
uncorrectable errors accumulate in the three-qubit code is proportional to 1/k in the Markovian case, and to 1/k^ 
in the non-Markovian case. These differences have the same cause, since the equilibrium weight in the error space is 
closely related to the rate of uncorrectable error accumulation. 

The quadratic difference in the error weight between the Markovian and non-Markovian cases can be attributed 
to the existence of a Zeno regime in the non-Markovian case. Regardless of the correlations between the density 
matrix inside the code space and the environment, if the lowest-order errors are correctable by the code, there exists 
a Zeno regime in the evolution of the code-space fidelity. The effective reduction of the error rate with the rate of 
error correction for non-Markovian error models depends crucially on the assumption that the time resolution of the 
continuous error correction is much shorter than the Zeno time scale of the evolution without error correction. This 
suggests that decreasing the time for a single (infinitesimal) error-correcting operation can lead to an increase in the 
performance of the scheme, even if the average error-correction rate goes down. 

While in this paper we have only considered codes for the correction of single-qubit errors, our results can be 
extended to other types of codes and errors as well. As long as the error process only produces errors correctable by 
the code to lowest order, an argument analogous to the one given here shows that a Zeno regime will exist, which leads 
to an enhancement in the error-correction performance. Unfortunately, it is very difficult to describe the evolution of 
a system with a continuous correction protocol, based on a general error-correction code and subject to general non- 
Markovian interactions with the environment. This is especially true if one must include the evolution of a complicated 
environment in the description, as would be necessary in general. A more practical step in this direction might be 
to find an effective description for the evolution of the reduced density matrix of the system subject to decoherence 
plus error correction, using projection techniques like the Nakajima-Zwanzig or the TCL master equations. Since 
one is usually interested in the evolution during initial times before the codeword fidelity decreases significantly, a 
perturbation approach could be useful. This is a subject for further research. 
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